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Abstract. Novel results for the self-consistent single-particle spectral function 
and self-energy are presented for non-degenerate one-component Coulomb systems 
at various densities and temperatures. The GVl^'^'-method for the dynamical 
self-energy is used to include many-particle correlations beyond the quasi-particle 
approximation. The self-energy is analysed over a broad range of densities and 
temperatures (n = lO^'^cm-^ - 10^'' cm'^, T = 10^ eV/fce - 10*eV/fcB). The 
spectral function shows a systematic behaviour, which is determined by collective 
plasma modes at small wavenumbers and converges towards a quasi-particle 
resonance at higher wavenumbers. In the low density limit, the numerical results 
comply with an analytic scaling law that is presented for the first time. It predicts 
a power-law behaviour of the imaginary part of the self-energy, ImE oc — n^/**. 
This resolves a long time problem of the quasi-particle approximation which yields 
a finite self-energy at vanishing density. 
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1. Introduction 

Strongly correlated Coulomb plasmas, fomid e.g. in planetary interiors [HE], fusion 
plasmas [3], and plasmas excited by lasers or ion beams [4J, are characterized by 
a high degree of spatial and temporal correlations, which lead to the emergence of 
phenomena like collective plasma modes, dynamical screening of the interparticle 
interaction potential, and dissolution of bound states. In particular, laser excited 
plasmas cover a broad range of densities and plasma temperatures. Values range from 
typical condensed matter conditions to hot, weakly coupled plasmas. 

Theoretical approaches to the physical properties of such systems have to deal 
with a great complexity. A particular challenge is the formulation of a coherent theory, 
which is valid over a wide range of densities (n) and temperatures (T) , thereby allowing 
to describe matter in various states, e.g. a solid-state target, being transferred into a 
plasma by interaction with high-power lasers and its subsequent relaxation [5j . Many- 
particle perturbation theory [6] presents a general approach to many-body systems like 
condensed matter partially and fully ionized plasmas and nuclear matter, to 
mention only a few. Also for non-abelian systems, such as the quark-gluon plasma [9] , 
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there exist similar approaches to the one described here for Coulomb systems, e.g. the 
concept of Schwinger-Dyson equations, see the review article [TD] . The thermodynamic 
properties as well as the response to external perturbations of these systems in various 
situations can be studied systematically [IX,. 

The central quantity within the many-body theoretical approach is the single- 
particle spectral function A{p,Lu). It represents a physical observable which can be 
measured via angular resolved photoemission spectroscopy (ARPES) [12l [131 El- 
Starting from the spectral function, a number of interesting questions related to 
the physics of many-particle systems can be addressed. The equation of state 
[15], transport cross-sections [T6j (e.g. electrical conductivity, thermal conductivity, 
and stopping power [17]), and optical properties [18] (emission and absorption of 
electromagnetic radiation) become accessible. 

In this work, the focus is on the spectral function of plasmas. As an example, 
a one-component electron plasma is considered which is charge compensated by a 
homogeneously distributed background of positively charged ions (jellium model). 
The plasma is characterized by the degeneracy parameters 9 and the plasma coupling 
parameter T which are defined as 

?i2(37r2n)2/3 ' AneokBT \ 3 J ' ^ ' 

Here, the electron mass m was introduced, /cb is the Boltzmann constant. In this 
work, we consider only non-degenerate systems, S> 1, i.e. the thermal energy k^T 
is large compared to the Fermi energy Ep = H'^{3TT^n)^/'^ /2m. 

The calculation of the spectral function becomes challenging in the regime of 
strong coupling, i.e. when the plasma coupling parameter becomes comparable 
or larger than unity. The coupling parameter measures the ratio of the Coulomb 
interaction energy of two particles at a mean distance to their thermal energy fceT. 
At r > 1, particle collisions become frequent, involving transfer of both momentum 
and energy. The interparticle potential is screened due to the presence of nearby third 
particles. These correlations significantly modify the plasma observables and have to 
be accounted for in the calculation of the spectral function. This is accomplished via 
the single-particle self-energy E(p, w), which is a complex function of both wavevector 
p and frequency oj, leading to a structured spectral function. Though, the main 
task of many-particle theory, applied to strongly coupled systems, is to calculate the 
self-energy in a suitable approximation. 

The simplest approximation, often found in the literature on Coulomb systems, 
is the mean- field or Hartree-Fock approximation [8]. One obtains a frequency 
independent self-energy which induces a shift in the spectral function's pole, the so- 
called Hartree-Fock or quasi-particle shift. For dilute plasmas, this correctly describes 
the lowering of the chemical potential due to the averaged field of the plasma particles. 
Also, the shift of the ionization energy for bound states is obtained [191 [20]. However, 
in dense systems, the mean-field approximation breaks down since the dynamical 
screening and collective excitations cannot be accounted for. One has to go beyond 
the quasi-particle picture. 

A particularly successful approximation for the self-energy, including these 
dynamical effects, is the so-called GW^-approximation [3T1 [52]. Correlations are 
accounted for via the dynamically screened interaction potential W{q,Lu), rather 
than via the bare Coulomb interaction. The GW approximation knows a long 
history of applications in the field of condensed matter theory. Examples are the 
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calculation of single-particle spectra in the homogeneous electron gas |23[ [Ml [5S], 
bandgaps in semiconductors effective masses of metal electrons [57], optical 

and electronic properties of insulators [28j . electronic structure of superconductors 
[25] . but also atomic and molecular systems [301 El In particular, GW self- 

energy corrections systematically improve band-gap calculations performed by means 
of density functional theory [33l |34l [35] . 

Recently, the GW approximation has been applied also to dense plasmas. 
Whereas Fehr et al. [SfT performed lowest order (one-loop) self-energy corrections 
to the equation of state, Wierling et al. [37] carried out pioneering self-consistent 
calculations of the electron self-energy in the solar-core plasma. An asymmetrically 
broadened, otherwise featureless spectral function was obtained. In this work, the GW 
self-energy and the corresponding spectral function is investigated for non-degenerate, 
one-component electron plasmas. Only unbound electrons are considered, bound state 
contributions can be accounted for via T-matrix calculations, as done in [35] . The self- 
energy is evaluated for a broad range of densities and temperatures, going from ideal, 
weakly coupled plasmas (F <C 1) to the strong coupling regime F > 1. As a novel 
contribution to the field, an analytic scaling law for the GW self-energy at low densities 
is derived which accurately describes the numerical data in this limit. This expression 
can be combined with corresponding formulae that are valid in the degenerate case, 
when k^T ^ E-p, to construct a fit formula for the self-energy which then covers a 
large portion of the density-temperature plane. 

Formerly, analytic expressions for the self-energy have been derived that base 
on the quasi-particle approximation [39]. In particular, the completely degenerate 
electron gas at T = was considered, using the plasmon-pole approximation |40j . 
and also weakly coupled (F <C 1), classical plasmas {Ep <C k^T), using the Born 
approximation for the self-energy [8]. The latter result exhibits several problems: The 
imaginary part of the quasi-particle self-energy is independent of density and carries a 
prefactor cc 1/h. Thus, there is an unphysical finite damping of single-particle states 
even in the vacuum and the classical limit /i ^ is not defined. On the other hand, 
from physical arguments, one expects that the self-energy vanishes at zero density 
and that it is a purely classical expression {h — 0), when 6^1. This problem has 
remained unresolved up to now. The real part of the quasi-particle self-energy is well 
behaved, i.e. it vanishes at zero density and is purely classical. 

The new analytic expression for the self-energy presented in this paper is derived 
without the quasi-particle approximation, i.e. it is a non-perturbative result. It is 
shown that only this non-perturbative treatment leads to an expression that is classical 
for both the real and the imaginary part and vanishes exactly in the vacuum limit 

The work is organized as follows: After a brief recapitulation of the single-particle 
spectral function and the GW-method in section [21 numerical results for the self- 
consistent spectral function and self-energy will be discussed in section [3] Section 3] 
contains the derivation of the non-perturbative scaling law and comparison to the 
numerical results. In section [5] it will be analyzed why the quasi-particle picture is 
incapable to give a physically consistent result for the imaginary part of the self-energy. 
Conclusions will be drawn in section [6] The appendix contains detailed calculations 
that are only summarized in the main part of the paper. 
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2. Spectral function and self-energy 

The derivation of the GW^-approximation involves some lengthy manipulations. In 
this section, only the most relevant formulae are given, while [Appendix A| contains 
the detailed steps. 

Central to the description of electronic properties in a many-body system, 
which is in thermodynamic equilibrium, is the thermodynamic electron single particle 
Green function G'(p,Zy), defined at the discrete Matsubara frequencies z,j = {2v + 
l)7ri k^T/h, 1/ = 0, ±1, ±2, .... It is related to the single-particle self-energy S(p, z,^) 
via Dyson's equation 

G(p, z,) = G(") (p, z,) + G(") (p, z,)E(p, z,)G(p, z,) 

G(")"'(p,z.)-E(p,z.)]"' , (2) 

with the free Green function G'-*'-' (p, z^) = [hz,^ — Sp] ^. Also, the single-particle 
energy Sp — It?p^ /2m — /i is introduced, fi is the electron chemical potential. G(p, z^) 
contains the thermodynamic properties of a single particle coupled to a thermal bath 
at a given temperature T. For example, the momentum distribution function is easily 
obtained by summation of the Green function over all Matsubara frequencies. 



n(p) = fcBr^G(p,z.) . (3) 

Instead of the complex Matsubara Green function, it is more convenient to operate 
on the real valued spectral function A(p,uj), defined on the real frequency axis. It 
carries the same information as the Green function and is defined via the spectral 
representation of the latter, 

G P,2. = / 7^ ^ ■ (4) 

Here, w is a real valued frequency. This relation can be resolved for A(p, w), 
A{p,io)= - lim 2lmGip,Lu + iS) (5) 

-2ImS(p,c^ + i^) 

= hm 5 5- , (6) 

•5^0+ [huj - £p - Re E(p, uj)Y + [Im S(p, uj + iS)Y 

i.e. the spectral function is obtained after analytic continuation of the Green function 
from the Matsubara frequencies to arbitrary complex frequencies as the imaginary 
part of G(p,a; + iS), when S approaches zero from positive values. In this way, the 
sign of the imaginary part of the self-energy is fixed, i.e. ImE(p(jj) < for (5 > 0. The 
real part of the self-energy behaves unambiguous for 5 = 0. 

The spectral function usually exhibits several resonances, including a central peak, 
located at the quasi-particle energy Ep, i.e. the solution of the quasi-particle dispersion 

= ep + ReS(p,£;p/;i) , (7) 

accompanied by symmetrically distributed satellites which are attributed to collective 
modes in the many-particle system (40j . The width of the resonances in the frequency 
domain is commonly identified with the inverse life-time of these excitations. 

Let us first look at the lowest order approximation to the self-energy, the Hartree- 
Fock term I]^^(p). It is given by the convolution of a non- interacting Green function 
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with the unscreened Coulomb potential V{q) = e^/eoq^Qo (Oq is a normalization 
volume), 

EHF(p^ = -kBTj] G(°) (p - q, - LO^) V{q) (8) 



q 



[l-nF(ep_q)]n9)^S"F(p), (9) 



with the Fermi distribution function n-pifuS) = [exp(?ia;/fcBT) + 1] . In the 
first line, summation takes place over the Bosonic Matsubara frequencies = 
27ri/x A;bT//i , /i = 0, ±1,±2,.... The first term (Hartree term) diverges, 

but it is exactly compensated by the same term from the positive charge background. 
The second term (Fock term or exchange term) gives a finite contribution. Closed 
expressions can be given in the case of non-degenerate plasmas [51 [32] and completely 
degenerate Fermi gases [7j. One finds E^^(p) cx n in the high temperature limit 
(fceT" 3> -Bp) and E^^(p) cx n^/'^ in the quantum degenerate case k^T <C i?F, see [5] 
for details. Thus, the Hartree-Fock self-energy fulfills the physical constraint to vanish 
at zero density. 

The Hartree-Fock term is a real function of momentum, only. The corresponding 
spectral function is shifted from the free particle dispersion, 

^HF(p,c^) = 27r^(ep + EHP(p)-;iw) . (10) 

No imaginary part of the self-energy appears in this approximation, i.e. the life-time 
of the Hartree-Fock quasi-particles is infinite. This is consequence of the mean-field 
approximation, where no fluctuations of the electric field, i.e. no dynamics of the 
surrounding plasma particles are taken into account. Recently, also the second order 
exchange contribution to the self-energy has been obtained in closed form [4ll [42] , 
see also [43|. However, this term and all higher order terms, involving only the 
bare Coulomb potential, do not lead to a finite particle life-time, only a shift of the 
dispersion relation is obtained. 

To describe the situation in a dense and strongly correlated system, where the 
single particle states are spectrally broadened, i.e. they acquire a finite life-time, one 
has to go beyond the quasi-particle approximation, and take into account the screening 
of the interaction. The GVF-approximation, can be regarded as the generalization of 
the Hartree-Fock theory to dynamically screened interactions. It was introduced by 
Hedin [53] for the homogeneous electron gas, and is defined as 

E(p, z,) = -k^T J2 G(P - q, z, - LU^) W{q, u;^) . (11) 
q.i^M 

W{q, z) is the dynamically screened interaction. Note that the GW approximation is a 
self-consistent ansatz, since the self-energy appears on the l.h.s. as well as in the Green 
function on the r.h.s. of (|lip . Also, the screened interaction W{(\,ljj^) is a functional 
of the Green function via the dielectric function e(q, o;^), i.e. the polarization fmiction 
n(q, uJ^j): 

In GW-approximation, n(q, uj^) is given by the inner product of two Green functions, 
H(q,a;^) — —k^TY^^ ^^G{q^ -\- \i, z^ -\- uo^j) G{Tp, z^). The "double" self-consistency 
implied in this ansatz makes the GH^-approximation is complicated and a numerically 



Spectral Function for Non-Degenerate Coulomb Systems 



6 



demanding problem. On the other hand, the fuh GW^-approximation suffers from 
deficiencies due to the neglect of vertex-corrections [H], such as violation of the /- 
sum rule 1451. This problem can be avoided by keeping the dynamically screened 
interaction on the level of the random phase approximation (RPA) [46j . defined by 
the RPA polarization function, 

nRPA(q,c^^) = -A:BTVG(°)(p-q,z,-c^^)G(°)(p,z,) , (13) 



nRPA(q,c. + i<5) = -V . (u) 



'^F(ek+q/2) - '^F(ek-q/2) 



k 



h{uj + i6) + ek-q/2 - £k+q/2 



The use of the RPA polarization function leads to the so-called GW^^°^ approximation 
for the self-energy. It has been shown to give more accurate quasi-particle energies 
[25] than the full GW-approximation. Additionally, it is known that higher 
order corrections beyond GW^°\ such as vertex-corrections and corrections in the 
polarization function beyond RPA, partially compensate. Therefore, ignoring them 
altogether is expected to give better results than accounting for one or the other [22] . 
The /-sum rule is fulfilled. Further sum rules, e.g. for the moments of the spectral 
function can be derived [24] which are useful to control the numerical treatment of 
the integral equations to solve. 

The inverse dielectric function epp^(q, cj) describes the propagation of 
electromagnetic waves in the plasma. As a main feature, it contains the longitudinal 
plasma oscillations or plasmons. These resonances show up as peaks in the inverse 
dielectric function, located at the roots of the plasmon dispersion Re eRpA(q, t^^) = 0. 
For non-degenerate systems, as considered here, the plasmon dispersion can be 
expanded in powers of the wavenumber q, and one finds the Gross-Bohm relation [47] 
^resil) = ^pi(l + l'^/^'^) + {fi^q^ f^m)"^ for the plasmon resonance frequency uiresiq)- 
Here, the plasma frequency ujpi and the inverse Debye screening length k 

-,1/2 

, (15) 



21 1/2 



eom 



have been introduced. A detailed discussion of the plasmon resonance in dense plasmas 
can be found in [48] . For the present discussion, it is important to keep in mind that 
the collective plasma excitations are accounted for via the inverse dielectric function 
in RPA. This is the main advantage of the GVF^°^-approximation compared to the 
mean-field or Hartree-Fock approximation. Depending on the choice of parameters 
like density and temperature, these plasmon resonances determine the shape of the 
self-energy as a function of the frequency and thereby also the spectral function, where 
satellites besides the quasi-particle peak indicate coupled electron-plasmon modes, 
often referred to as plasmarons [1^ . 

It should be noted at this point that contributions from bound states to the self- 
energy are not accounted for in this work. The description is limited to fully ionized 
plasmas. Bound state contributions can be included using the concept of the T-matrix, 
see e.g. the work by Schmielau et al. [35]. 

Using the spectral representation ^ and the screened interaction p^ . the 
following equation for the imaginary part of the self-energy in GVF'^*'-'-approximation 
is obtained after summation over the Bosonic Matsubara frequencies tj^. 



Pi /"OO J / 

Im S(p, + i^) = / ^V{q)A{p - q, - cj') 

nF[hu;) q 7-cx) 27r 
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xlmej^p^{q,ij') nQ{hw') np{huj — huj') , (16) 

with the Bose-Einstein distribution function n-Q{huj) = [cxp{hijj / k^T) — 1] ^. The 
real part of the self-energy is obtained by means of Hilbert transform as 

ReE(p,.) = ESr(p) +7^ r ^ ^^^^^ . (17) 

V denotes the Cauchy principal value integration, S]P^(p) is the Hartree-Fock self- 
energy of the interacting system, 

/OO 1 
^^A(p-q,c.)np(M^^('Z) ■ (18) 

Finally, to close the set of equations, the chemical potential has to be fixed by inversion 
of the density relation 

n(M,T) = 2— 2^A(p,c^)nF(M . (19) 

The factor 2 in front of the r.h.s. stems from the summation over the spin components. 
Together with Dyson's equation ([6]), (fT6|) - p9|) constitute a system of non-linear 
integral equations for the self-energy. 

Besides the normalization of the spectral function 



(20) 



similar sum-rules can be derived also for higher moments of the spectral function 
|24| . These are independent of the concrete approximation used for the self-energy. 
In second order, one obtains an equation relating the first moment of the spectral 
function to the interacting Hartree-Fock self-energy (|18p . 

^u;Aip,cu)=e^ + j:fJip) . (21) 
, 27r 

Similarly, the second moment is related to the Hartree-Fock energy and the frequency 
integrated imaginary part of the self-energy, which is itself a conserved quantity, at 
least within the GW'^°^ approximation, see (j23p below, 

(22) 

For the GW^^"^ self-energy, Holm and von-Barth have found the following identity, 
relating the integrals over the imaginary part of the self-energy to the totally integrated 
response function. 



n — Iml](p,^ + i5) = ft^ / —V{q)lmepl,f^{(i,L 



(23) 

q 

In the next section, results for the self-energy will be presented that are obtained 
via numerical solution of (|16p . The sum rules given above are used to check the 
accuracy of the numerical results. 
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3. Numerical results 

The GW^'^'-approximation is evaluated numerically for various sets of plasma 
parameters in the following. A typical example of a weakly coupled (F = 0.07), 
moderately degenerate {9 = 2.2) plasma is the plasma at the solar core, with 
temperatures of T ~ lOORy/fce — 1360eV/fcB and electron densities of n ~ 
7 X 10^^ cm~^ [49]. The solar core plasma has been investigated using the GW^*^"^- 
method in a number of previous publications, see [IHl [SZl HO] ■ Here, most attention is 
paid to a systematic analysis of the single-particle spectral function and the self-energy 
over a broad range of densities and temperatures, however, sticking to non-degenerate 
plasmas and neglecting bound states. We therefore start with a plasma temperature 
that equals the solar core temperature and a density that is 10% of the solar core 
electron density. Later, higher and lower temperatures will be considered as well, i.e. 
kBT = 10 Ry and /cbT = 1000 Ry. Note that is always chosen large against 
typical binding energies of atoms which are usually of the order of several Ry. Thus, 
bound states can be neglected. 

The numerical solution of equation (jl6[) is performed by means of an iterative 
algorithm, starting from a suitable initialization of the spectral function. Typically, 
the algorithm converges after 5-10 iterations. The threefold integral (fT6|) is evaluated 
on a two dimensional grid with roughly 100 nodes in the frequency coordinate and 
10-20 nodes in the momentum coordinate. The angular integral is performed first, 
followed by the frequency integration and the integration over the modulus of the 
wavenumber q. The result is checked for consistency in each iteration using the sum- 
rules pm23|) . Further details concerning the numerical implementation are provided 
in [SI]. 

Figure d] shows the numerical result for the self-energy (imaginary and real part), 
dispersion relation ftw-f — fi,^p^/2m — Re I](p, cj), and the spectral function for plasma 
parameters chosen as rt = 7 x lO^"' cm^'^ for the plasma density, i.e. 10% of the solar 
core density, and T — 100 Ry = 1360 eV/ /cb for the plasma temperature. The chemical 
potential is ^ — — 377Ry = — 5133eV. The momentum was fixed at hp = 0. 

The spectral function at the chosen parameters is a broadened resonance with 
two satellites appearing at about hio + fJ, = ±5 Ry which is slightly below the plasma 
frequency hiOpi = 7.2 Ry at the chosen conditions. The plasma frequency is indicated 
by the dashed vertical lines. As already mentioned, these satellites are often referred 
to as plasmarons, i.e. a coupled mode between the single particle resonance and the 
collective plasma oscillation [40]. The imaginary part of the self-energy (top graph) 
is peaked at the free dispersion huj + ij, = Q and this peak leads to the small dip in 
the spectral function between the satellites, sec The real part of the self-energy 
(second graph from top) is a rather smooth function, leading to only small variations 
in the dispersion (3rd graph from top). 

Next, the dependence of the spectral function on the wavenumber p is analyzed. 
In figure 121 the spectral function A{p,u!) is shown for five different wavenumbers, i.e. 
p = 0, 5ag^, lOflg^, ISflg^, and 20 a^^, ob = Aneoti^ /me^ is the Bohr radius. The 
density and temperature are the same as before, n = 7 x 10^''cm~'^, /cbT = 100 Ry. 
At increased wavenumber {p > 5a^^), enhanced complexity of the spectral function 
is observed. The plasmaron peaks, which at hp = appear as small shoulders 
in the otherwise broad central resonance, are better defined. The central quasi- 
particle peak itself becomes narrower and the plasmaron peaks separate. At the 
highest momenta considered {hp > 15 a^^), the plasmarons themselves are damped 
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energy ti[o+|j, [Ry] 



Figure 1. Self-energy (in units of the Rydberg energy, 1 Ry = 13.6 eV), 
dispersion relation (in units of Ry), and spectral function (in units of 1/Ry) 
for plasma density n = 7 X 10^* cm""^ (10% of the solar core density) and 
temperature T = lOORy/ZcB = 1360eV/fcB. The spectral function contains two 
weakly pronounced plasmaron satellites, appearing at slightly smaller energies 
than the plasma frequency (dashed vertical lines). The chemical potential is 
H = -377 Ry = -5133 eV. 

out, and a single, narrow resonance forms, focated near the single particle energy 
hw = Ep = h^p^ /2m — ii,i.e. the quasi-particle picture is restored. Some of these 
features, especially the plasmaron satellites are already known from literature j40| . 

Now that the general characteristics of the spectral function have been discussed, 
the central concern of this paper can be worked out, i.e. the analysis of the dependence 
of the self-energy and the spectral function on the plasma parameters density and 
temperature. In figure [31 the spectral function at p = is shown for five different 
densities between n — 7 x 10^^ cm~^ (solar core conditions) and 0.01% of the solar 
core density. The temperature is kept constant at T = lOORy/fce = 1360eV/fcB- 
The spectral function drastically changes with varied density. In the case of the 
highest density considered, a narrow quasi-particle peak accompanied by two separate 
plasmaron satellites (indicated by arrows) is observed. The quasi-particle peak is 
notably shifted from the free dispersion Eq = fM, due to the real part of the self-energy. 
Going to lower densities, the plasmaron satellites merge into the central peak, as can 
be seen in the case of the spectral function for n = 7 x 10^^ cm^^ and also the quasi- 
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Figure 2. Spectral function for plasma density n = 7 X 10^* cm~'^ and 
temperature k^T = 1360 eV (solar core temperature) as a function of momentum 
and density. The black line on the bottom represents the free dispersion relation 
hui = £p = T?'p'^ jlm — ji. At the present parameters the chemical potential is 
^ = -377 Ry. 



Pi 



o 
II 



n=7.0xl024 cm'^ 
n=7.0xl022 cm'^ 
n=7.0xl022 cm'^ 
n=7.0xl02j cm'^ 
n=7.0xl0 cm 




-20 20 

energy hco+iJ, [Ry] 



Figure 3. Spectral functions at p = for different plasma densities, ranging from 
the density at the solar core, n = 7 X Vfi^ cm~^ (solid curve) to 0.01% of the 
solar core (dash-dot-dotted curve). The plasma temperature is T = lOORy/fee = 
1360eV/fcB for all five curves. 



particle shift is reduced. Finally, at the lowest densities considered, n ~ 7 x 10^^ cm^"^ 
and 7 x lO^^cm^'^, a single, narrow quasi-particle resonance is obtained which is 
centered around the bee dispersion. The width decreases with the density which is 
the expected behaviour in the low density limit. 

In order to study the dependence of the self-energy on density and temperature 
in more detail, the effective quasi-particle self-energy I]{p, Ep/h) as a function of 
the density at various temperatures is considered. This quantity gives the shift and 
width of the central peak in the spectral function ^(p,^;), i.e. when lu is close to 
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Figure 4. Effective quasi-particle damping width at p = 0, —lTaT,{0, Eq/K) 
normalized to the thermal energy as a function of the plasma density. The arrows 
indicate for each temperature the density at which the degeneracy parameter 9 
takes the values 6 = 10 and 9 = 1. 



the quasi-particle frequency Ep/h, see (O. The resuhs for the imaginary part of the 
effective quasi-particle self-energy at p = as a function of the plasma density are 
shown in figure 21 Three different temperatures have been assumed, T = 10, 100, and 
lOOORy/fee- Towards low densities, a systematic decrease of — ImE(0, Eq/K) with the 
density is observed which is also known from the literature |37j . The asymptotes to the 
low density behaviour, shown as thin dotted lines, indicate that Im I](p, iJp/?i) scales 
proportional to —n^^^. This behaviour will be analyzed in more detail in section [H 
where an analytic solution for the GW self-energy is derived that exhibits the same 
v}/'^ proportionality. 

At higher densities, the power law behaviour terminates and the self-energy starts 
to decrease. This can be understood by looking again at figure [3l Here, it was shown 
that at increased density, the plasmaron satellites separate from the central quasi- 
particle peak, i.e. spectral weight is shifted to the satellites and the central peak 
narrows. The calculations have only been performed for non-degenerate systems, i.e. 
for densities, where the degeneracy parameter 9 = k^T / E-p is still large compared 
to unity. The extension to degenerate systems is straightforward and is covered in 
another paper [51], but will not be treated in this work. 

The real part of the self-energy (effective quasi-particle shift), at the densities and 
temperatures considered here, was found to follow exactly the Hartree-Fock behaviour, 
i.e. ReS(0, Eq/H) = ~h?K^ /2m (x —n [8]; k is the inverse Debye screening length, see 
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4. Analytic solution for the GW'^^^ self-energy in Born approximation: 
Classical limit 

4.I. Derivation of the analytic solution 

As discussed, the spectral function in the low density limit is lacking any plasmaron 
resonances, only a broadened quasi-particle peak appears, see figure |3l In order 
to understand this behaviour, the GW^^*^^ -equation is reconsidered applying a 
sequence of approximations as described in the following. In this way, an analytic 
solution is found that is valid at low coupling parameters. 

It will be shown that the observed scaling is obtained correctly, if the imaginary 
part of the self-energy is kept finite also on the r.h.s. of the self-energy integral 
equation (fT6|) . It therefore represents a generically non-perturbative result. Details of 
the calculations can be found in [Appendix B| 

Since collective excitations do not show up in the self-energy and the spectral 
function at low densities, it is obvious to neglect these features already in the screened 
interaction. Formally, this is achieved by replacing the complete inverse dielectric 
function by the Born approximation. 



Ime '(q,u;) ~ — , , 



e^q,0)P 



(24) 



For the static dielectric function appearing in the denominator, we use the Debye 
expression eD(q, 0) = 1-1- K^/q'^, with the inverse Debye screening length. In other 
words, instead of the interaction via a dynamically screened potential, electron- 
electron collisions via a statically screened potential are considered using the Born 
approximation. Then, l|16p turns into 

Imi;(p, w + iS) = 



2mkBTe'^K^ f\ r, f°° . , f°° dgg / mtj'^ \ f hoj 

dcose/ doj 2 *5xp — — — exp 



^3 4^eoy_i J-oo Jo (q^ + K^y V '^q^ksTj " VSfceT 

Ini E(p — q, w + i(5 — ui') 
[tiw — fku' — £p-q — Re I](p — q, — uj')]'^ + [Im I](p — q,uj + iS — uj')]^ 



(25) 



Note that the dielectric function is taken in the classical limit, i.e. the Fermi-Dirac 
distribution is replaced by the Maxwell distribution, leading to the exponentials in 
the first line of 

Due to the statically screened Coulomb potential, important contributions to the 
g-integral stem from values q ^ k. Therefore, we neglect the shift of momentum in the 
self-energy on the r.h.s. of equation ((25)) . i.e. we write S(p--q, oj — oj') ~ T,{p,uj — u!'). 
To justify this approximation, we show the numerical solution for the imaginary part 
of the self-energy (P5|) in figure [5] (dashed curve). The solid curve corresponds to the 
solution that is obtained by neglecting the momentum shift in the argument of the 
self-energy on the r.h.s. of (P5|) . As can be seen, this approximation docs not modify 
the result significantly. In fact, the small deviations, which are only observable around 
Hlu -|- /i ~ 0, are already in the order of the numerical accuracy. 

Subsequently, the remaining terms in (|25p are expanded in powers of q/k, as 
described in detail in the appendix. Finally, the threefold integral can be performed 
and the equation 

2 

[ImE(p,w + i(5)]^ + r?iV/2TO-Ai - ^ + RcS(p,tj)]^ = fceT-^ (26) 

47reo 
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Figure 5. Imaginary part of the self-energy at momentum hp = for plasma 
parameters n = 7 X 10^^ cm""^ and T = lOORy/fce- The self-consistent Born 
approximation (finite q, dashed curve) is compared to the calculation where the 
momentum shift hq is neglected in the self-energy on the r.h.s. of the self-energy 
equation (solid curve). 



is obtained. The l.h.s. is just the denominator of the spectral function, c.f. ([6|). 
Together with the spectral representation of the Green function ^ , we then find the 
equation 



[ftz - ;iV/2m + /i - S(p, z)] ^ = 



which, in the limit z = uj + i6, S 0^ has the solution 
E(p, uj + i5) = 



S(p,z) , 



—sign{huj — h^p^/2m + /i) 

The signum function, 

sign(w) = 



huj + i5 — hrp^ /2m + ^ 



-,1/2 



47reo 



fceT 



1 ^ w > 
1 ^ u <0 



(27) 



(28) 



(29) 



ensures the correct sign of the imaginary part of the self-energy, i.e. Im I](p, a;+i(5) < 
for (5 > 0. 



4-2. Comparison to the numerical solution 

The imaginary part of is plotted in figure [5] for T = lOORy/fce and n = 
7 X 10^^ cm"'^, i.e. for the smallest density considered in figure ([3]). The analytic 
formula is compared to the full numerical solution for two different wavenumbers, 
p = {a) andp = I/ob (b). The dotted vertical line indicates the position of the quasi- 
particle dispersion Ep. In the first case, both numerical and analytic calculation agree 
reasonably well, albeit the analytic solution lies systematically above the numerical 
data. However, the overall deviation is smaller than 7%. In the second case {p = I/ob), 
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the upshifted plasmon peak, present in the numerical result, is not reproduced by the 
analytic formula. Thus, the analytic formula is applicable only for small momenta, 
while at higher momenta, the dynamical features of the interaction become important. 
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Figure 6. Imaginary part of the self-energy for plasma density n = 7x 10^^ cm~'^ 
and temperature T = lOORy/fcB- Results for p = (a) and for p = l/ae (b) are 
shown. The self-consistent GVyC'-calculation (solid curve) is compared to the 
analytic formula II28I I given as dashed curve. The dotted vertical line indicates 
the quasi-particle dispersion fkj = Ep. 

On the other hand, the analytic formula is very useful to initialize the numerical 
algorithm. This is analyzed in figure [T] Here, the spectral function, that is obtained 
in the first iteration of the algorithm, was computed in two different ways for the 
same parameters as above, n = 7 x 10^^ cm""^ and /cbT ~ 100 Ry. The dashed curve 
gives the first iteration starting from the analytic formula (|28|) for the self-energy, the 
dotted curve is the same calculation but starting from a narrow Gaussian spectral 
function with a width of 0.3 Ry (FWHM). In plot (a) the wavenumber is p = 0, 
while in (b), p = l/ae was chosen. For p — 0, the analytic ansatz leads to a good 
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resemblance with the converged result (solid curve). The converged result is taken 
here as the 20. iteration starting from the Gaussian ansatz. The Gaussian ansatz, 
iterated once, results in a two-peak structure which is far from the converged spectral 
function. Also at p = I/ob, starting from the analytic ansatz gives a much better 
overall correspondence than the calculation starting from a Gaussian spectral function, 
although subtle details like the plasmaron peak at fex; + ~ 4 Ry is not reproduced 
in the first iteration. 
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Figure 7. Spectral function for plasma density n = 7 x 10^^ cm~^ and 
temperature T = 100 Ry /fee- Plot (a) shows the spectral function at p = 0, in 
(b) p = l/ao was chosen. The first iteration starting from a sharp quasi-paxticle 
spectral function (dotted curve) is compared to the first iteration starting from 
the analytic expression for the self-energy (dashed curve). 

In order to perform a quantifiable comparison between both initializations and 
their impact on the convergence of the algorithm, we determine the mean squared 
deviation of the spectral function in a given iteration u from the converged result 
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5^ = AT-i J2iLi {A^"^ (0, uji) - ^(20) uj.,)) with N the number of points on the uj- 
grid of the spectral function, LUi the grid points. The resuh is shown in figure [8] for 
the Gaussian ansatz (marked +) and the initiahzation using the analytic self-energy 
(marked x). During the first four iterations, the mean squared deviation of the second 
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Figure 8. Mean square deviation of the spectral function at p = from the 
converged result as a function of the iteration number u. Red symbols correspond 
to results obtained when starting with a sharp quasi-particle spectral function, 
the green symbols result from starting with the analytic self-energy and leads 
to relatively small deviations already in the first iterations. High accuracy 
{S < 10~^) is obtained in both schemes only after about 10 iterations. Plasma 
parameters: n = 7 X lO^^ cm'^, T = lOORy/fce- 



method is by two orders of magnitude smaller than if using the Gaussian spectral 
function. While the mean squared deviation using the analytic self-energy becomes 
smaller than 10"^ already after 5 iterations, it takes 8 iterations for the Gaussian 
ansatz to get to this point. Also, it was found that a Gaussian with the width fixed 
at the imaginary part of the effective quasi-particle self-energy, does not improve the 
convergence, since the special form of the self-energy and the spectral function with 
a broad plateau and steep edges cannot be reproduced by such an ansatz and the 
analytic self-energy given in (^5)1 should be used instead. 



4-3. Analytic solution at the quasi-particle dispersion 



In the following, the analytic solution ([28)1 with the frequency fixed at the quasi- 
particle dispersion oj — Ep/h shall be considered in more detail. The numerical 
results for T,{p, Ep/h) at p — have already been discussed in section [3l see figured] 
Since the only dependence on frequency and wavenumber is given by the trivial term 
hu! — h^p^/2m -\- fi = hu! — Sp, the discussion may be restricted to the case p = and 
hu! = Eq = —ji. Note that due to ([^5]) Rel](p,£p) = 0, therefore Ep = Sp. Then, the 
imaginary part of ((28)) reads 



47ren 



47r n k-oT 



1/4 



(30) 
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which can also be given in terms of the plasma couphng parameter F, 

ImS(0,-M/?i) = - (ar^)^/'' /cbT . (31) 

This value is the damping width neglecting the influence of electron-plasmon 
interaction via the dynamically screened potential. In the following, it is referred 
to as the non-collective damping width. The non-collective damping width depends 
solely on the temperature and the classical coupling parameter F. Therefore, it is a 
purely classical result. 

We compare this result to the numerical solution for the GW^'^'' self-energy, 
presented in section [3) In figure [HI the effective quasi-particle self-energy at vanishing 













1 1 1 1 






1 y 


































10 Ry 










CO 
-i<! 


10-' 






- v= 

- v= 


100 Ry 
1000 Ry 




~+^- 


-0.5 


11=1 _ 


\ 














1 


















1 






o 
W 

S 

1 










1 


1 

71-0.5 








10-^ 






! 

Tl=0.5 


1 

11-1 

, , , , 1 






1 





t \ \ \ \ I \ 

-3 -2 -1 

10 10 10 



r 



Figure 9. Effective quasi-particle damping width at p = 0, Im S(0, -Eq/^) 
normalized to the thermal energy as a function of the plasma coupling parameter 
r. The solid black line marks the derived scaling law following II30I I. The vertical 
lines mark for each temperature the onset of plasmon excitation and corresponding 
decrease of quasi-particle damping, i.e. the parameter r) = ^Tjob > 1- 



momentum lml](0, i?o/^), in units of the plasma thermal energy, is shown as a 
function of F. The numerical results perfectly agree with the derived scaling law in 
the low density limit. As long as the density is small, such that the plasma frequency 
is below the non-collective damping width, i.e. 

ry^ , «1 (32) 

K e^kBT/ATTeo 

the damping is mainly non-collective, and a single, broadened resonance appears in 
the spectral function, c.f. figure [3l 

The parameter rj can also be expressed through the plasma coupling parameter 
F and the degeneracy parameter 9, or through the inverse Debye screening length k, 

1/6 

r-1/4 5,-1/2 _ Q gggg p-1/4 ^-1/2 ^ ^33-) 



V 



35/27r 



77 = ^/kob ■ (34) 

Thus, the derived scaling law is only valid for large 9, i.e. classical systems. This was 
already shown in figured When 9 approaches 1, quantum effects set in. For example. 
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collisions become less probable due to Pauli blocking which leads also to a decrease of 
the self-energy. 

Obviously, the Bohr radius sets the relevant length-scale that is to be compared 
to the inverse screening length k in order to estimate the importance of non-collective 
damping. Non-collective damping is the dominant mechanism, as long as the screening 
length is large compared to the Bohr radius. When the screening length becomes 
smaller than the Bohr radius, i.e. 77 > 1, which, due to ([32|) . is equivalent to 
having the energy of plasma oscillations larger than the non-collective damping width, 
the plasmaron satellites begin to separate from the broadened quasi-particle peak. 
Spectral weight is transferred from the wings of the central peak into the plasmaron 
satellites leading to a more defined quasi-particle resonance, i.e. a decreased damping 
of the central peak, see figure [3l Concluding, the analytic result is only a good 
approximation at low densities, when rj 1. At higher densities, dynamical screening 
becomes important, leading to satellites in the spectral function. 



5. Deficiencies of the quasi-particle approximation 

The non-collective damping width was introduced above in pop as the value of 
the imaginary part of the self-energy at vanishing momentum and frequency, p = 
0, oj + p,/h = 0. Of course, the same result is also obtained if this choice of variables 
was already made at the very beginning of the calculations leading to ((28l) . However, 
in the latter case, the manipulations can be performed in a different manner. At an 
intermediate step of the calculation, one can identify the reason why the quasi-particle 
damping Iml](p, Ep/h) as given in [HI [31], behaves unphysical in the low density and 
classical limits. Detailed calculations are given in [Appendix C[ while here only the 
most important steps are summarized. 

Setting p = and hus + p — in ([16|) . neglecting the momentum shift in the 
argument of the self-energy on the r.h.s. and replacing the dynamically screened 
potential by the statically screened Born approximation as before, we obtain 



Iml] 0,-MA = Re exp -z^ erfc iz , 35 

with 

h^q^ + i2mlm^{Q,-p/h) 

z = ^ . (36) 

2hq\/2m k^T 

Most contributions to the integral stem from small values of the wavevector, 
g < K. Therefore, we may neglect the real part of z and write z = 
i-\/2mImI](0, — n / K) / 2hq^J k^T . Using the expansion 

lim exp(x2) erfc(a;) ^ — i ^ 1_ „ + ©(x^^) (37) 

in lowest order only, the q-integral can be performed, resulting in 



Iml](0,-A*/^) = -V^^e^A^BTyiTfen , (38) 

which coincides with (j28p at p = and hw = —^l. 

From (|35p , one can also derive the quasi-particle approximation for the imaginary 
part of the self-energy: If the imaginary part of z, i.e. the self-energy, is neglected on 
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the r.h.s. (this is just the quasi-particle approximation), and furthermore the hmit 
g — > is considered, the expression 



is obtained. 

This coincides with the formula for the imaginary part of the quasi-particle self- 
energy as given in [51 p. 114, equation (4.164)]. There, the spectral function on the 
r.h.s. of the integral equation for E(p,a;) is replaced by an on-shell delta distribution 
(free particle spectral function), i.e. the self-energy is set to on the r.h.s.. The 
resulting integral is evaluated at the free particle dispersion fku = Sp. 

As result, one obtains the expression 



with iFi{a, (3; z) being the confluent hypergeometric function [52j. Note that in the 
given reference, instead of the imaginary part of the self-energy, the quasi-particle 
damping r{p,ep/h) = — 2Im I](p, ep/?i) is given. Also, the original formula differs 
from (|40p by a factor of 1/4. However, the formula given here was approved through 
private communication by W.-D. Kraeft. 

Obviously, 1^0)) is independent of density. The neglect of ImS(0, —fi/h) in the 
complex variable z leads to a different analytical structure of the equation. Therefore, 
the quasi-particle approximation has no chance to ever obtain the correct behaviour at 
low densities. Low densities, and therefore small inverse screening lengths n shift the 
supporter of the g-integral to small q, where contributions from Imz are important, 
whereas the real part of z vanishes at q — and leads to a result which is independent 
of K. 

In the same way, one can understand why the quasi-particle limit diverges when 
considering the classical limit ?i — s- 0. The imaginary part of z has h in the denominator 
which after the integration cancels the h in the prefactor in equation (j35p . No 
cancellation takes place, if the imaginary part is neglected, i.e. in the quasi-particle 
approximation. This leads to the divergence of the final result. 

6. Conclusion 

In this work, the single-particle self-energy of the one-component electron plasma was 
investigated. The spectral function was calculated self-consistently using the GW'-°^- 
approximation which allows the systematic treatment of dynamical correlations in the 
plasma. The spectral function contains at small momenta a broadened quasi-particle 
peak and two plasmaron satellites which, at low densities, merge into the central 
quasi-particle resonance. At increased momenta, for a given density and temperature, 
the spectral function converges to a single, sharp quasi-particle resonance. Special 
attention was paid to a systematic investigation of the self-energy and the spectral 
at different densities and temperatures. Here, only non-degenerate plasmas were 
considered, i.e. the temperature is large compared to the Fermi temperature. Also, 
bound states were neglected. 




(39) 




(40) 
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It was found that at low densities, the imaginary part of the on-shell self-energy, 
i.e. the inverse single particle lifetime, follows a universal scaling law ImS(p, Ep/H) cx 
— n^/^. For the first time, an analytic result for the on-shell single-particle self-energy 
was found that contains the correct low-density limit, i.e. a vanishing self-energy 
at n = 0. This is a major progress compared to the well-known quasi-particle 
approximation that yields a finite damping width even at zero density. The new 
on-shell single-particle damping width is — Iml](p, Ep/h) = (SF'^)^^* fesT. Since it is 
derived in Born approximation, i.e. no collective excitations contribute to the damping 
mechanisms, this quantity is called the non-collective damping width. By comparison 
of the numerical results to the new analytic formula, the parameter r] = ^/kob was 
identified to separate the regime of non-collective damping (77 ^ 1) from the regime, 
where the coupling between single particle states and collective excitations dominantly 
determine the single particle damping (77 3> 1) at small momenta. This analysis 
complements earlier work on the electron spectral function based on the plasmon-pole 
approximation. 

For 77 <C 1, the analytic formula (|28p is a good approximation for the self-energy. 
Furthermore, the use of the analytic formula for the self-energy as an initialization of 
the iterative algorithm leads to significantly faster convergence as compared to other 
methods, where a Gaussian ansatz is used as the initial spectral function. 

The non-collective damping is a purely classical result, no powers of h appear. 
This is fundamentally different from the quasi-particle approximation to the imaginary 
part of the self-energy which has no classical limit, i.e. the self-energy diverges in the 
limit /i — > 0. It could be shown that this problem, as well as the paradox of being 
density independent, stem from the inherently inconsistent treatment of the self-energy 
in the quasi-particle approximation. The long-time open question of the classical limit 
of the single-particle self-energy can now be regarded as settled. 

The results reported in this work are of paramount importance for many-particle 
theory and applications to dense plasmas. In particular, simple analytic expressions 
for the single-particle spectral function and self-energy in the classical and in the 
degenerate limit are needed to construct Pade-like interpolation formulae that cover 
the complete density-temperature plane. Such expressions would greatly simplify 
the calculation of equation of state, transport and optical properties of dense, high 
energy plasmas, solid state devices but also nuclear, hadronic, and partonic matter, 
and provide benchmarks for numerical approaches, i.e. simulation techniques. One 
part of this task, the analytic formula for non-degenerate dilute plasmas has been 
accomplished in this work. 
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Appendix A. Details on the GVF^^-'-approxiniation 

Throughout the appendix, the Rydberg system of units wiU be apphed to keep the 
formulae short and readable. In these units ?i = fee = 1, e'^ = 2, eg = l/47r, and 
m = 1/2. 

We start from the representation of the self-energy in terms of the full Green 
function G(p, z^), the dynamically screened potential W^(q, w^) and the vertex 
function r(p, p + q; z,y , + w^), given by the diagram 



W 



E(p,z^) = 




(A.l) 



p(o) G r 



In the GVF-approximation, the vertex is replaced by the bare vertex T^^^ — e, i.e. the 
charge of the considered particles, electrons in this case, 

W 




(A.2) 



p(o) G p(o) 



= -r^G(p-q,z,-c^^)W^(q,o;^) , (A.3) 

which is equation pip . The dynamically screened interaction is taken in the random 
phase approximation 

W-(°)(q,c.,) = (A.4) 

eRPA(q, w + i(5) = 1 -y(g)nRPA(q,w + i(5) (A. 5) 

nRPA(q,t^ + i(5) = - > — (A.6) 

^ UJ + ld+ £k-q/2 - 

Using the spectral representations of both the Green function ([4]) and the screened 
interaction in RPA, 

W^°\q,z)^Viq)(l+ r dc^Im^^yq^c^X ^^^^^ 



leads to 



doj A{p — q, w) 



q,i^M 

f^r- ^ime-^U^,.'y ^^^^ 



and after summation of the Bosonic Matsubara frequencies, 

—A{p-<i,u;") 

-oo 
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(A.9) 

is obtained. This expression contains the Hartree-Fock self-energy of the interacting 
system, 

/OC 1 
^^A(p-q,c.)nFMT/(g) , (A.IO) 

and the correlated self-energy 

^A(p-q,c.") 

^ r dc^- Imei^j>A(q.^O ["b(c^O + 1 - uAlo")] _ ^^^^^ 

J -00 I" - CJ' - CJ" 

For convenience, we skip the upper index "corr" in the following and only distinguish 
between the frequency dependent self-energy Yj{\),uj + iS) and the Hartree-Fock term 
E?^^(p), in the following. 

After analytic continuation Zi, — > z = uj + i6,S 0, the imaginary part of 
the correlated self-energy is evaluated using Dirac's identity lima^ol/C^; ± iS) = 

ImS(p,c^ + iO+) = — - V / -^V,MMp-^:^-^') 

xl-me^pj^{c[,uj')nB{uj')nF{uJ ~ Lj') , (A. 12) 

where the exact relation nB(^') + 1 — np{uj — lu') — —n-B{Lo')np{uj ~ uj')/n-p{Lo) was 
used. This equation is given as (|16p in the main text. 

Appendix B. Analytic self-energy for the classical one-component plasma 

In the high temperature limit k^T ^ Ep, we replace the Fermi Dirac distributions in 
the self-energy equation (fTB]) as well as in the dielectric function by the Maxwell- 
Boltzmann distribution, np{ek) — > f{k) = ^^exp(— et/T) with the thermal de- 
Broglie wavelength A — (47r/T)^/^. In this approximation, the polarization function 
takes the form [Ml, 



RenRPA(q, 



2qT 



to q 



2qVT 2Vt 



ImnRPA(q,w' 



^+,),fJi,3/2,-(^ + ^-'^ 



(B.l) 



Snq 



exp , , , 

2qVT 2VT, 

I t oj q \ 

exp 



2qVT 2VT 



(B.2) 



Spectral Function for Non-Degenerate Coulomb Systems 



23 



Then, the imaginary part of the self-energy writes 



Iml](p, w + iS) 



2k^ 
1 exp 



dcos6' 



|6(g,a;')P 



exp 



2T 



Im I](p — q, + ii5 — w') 



[a; — a;' — £p-q — Re I](p — q, a; — cj')]^ + Piii ^(P ^ ^ ^0]^ 

(B.4) 

Furthermore, the Born approximation is apphed, i.e. the dielectric function in 
the denominator is replaced by the Debye expression, £0(9,0) = 1 + n^jq^. 
Diagrammatically, the self-energy in this approximation is written as 




(B.5) 



G 



Vb denotes the Debye potential Vb(<z) — e^/eo((7^ -t- K^)rio. 

Since the main contribution to the q— integral stems from momenta q < k, we 
neglect the transfer wavenumber q in the argument of the self-energy on the r.h.s. and 
write 

2 2 j-l j-oa pea 

Iml](p,cj + i(5) = ^y^\/T / dcos6l dq dw' 

' J-l Jo J -00 



1 exp 



4(j2T 4T 



Im S(p, a; + i(5 — cj') 



[bj — u)' — ep-q — Re S(p, w — w')]^ + [Im S(p, w — t^')]^ 



(B.6) 



Furthermore, we neglect the term q^ /AT in the exponential which is small for high 
temperatures and for q < n, 



9 2 j-l j-oo noo 

ImS(p,cj + i5) = ^-^\/T / dcos6i / dq / duj' 

' J-l Jo J ~oo 



[g2 



exp 



UJ 



Aq^T 



exp 



LO 

2T 



Im S(p, w + i(5 — ijj') 



[uJ — Uj' — ffp-q — Re I](p, — 1.^')]^ + [I™ ^(Pj ^ ^ ^')]^ 



(B.7) 



This is equation was given in section 2] as ([^5)1 . 

Now, the integration over the angle can be performed as 



Im bj + i5) 



2ie 

7r3/2 



dcosfl 



dq 



duj' 



[q^ 
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w - w' - - + - Re w - w') _^ cog ^ f Im S(p, u + iS - tu')^ ^ 



2pq J \ 2pq 

ip + qf - uj - jjL ^-Re E(p, w - w') 



arctan 



arctan 



Im E(p, + i(5 — cij') 
{p — qY — (jj — ^ + ReS(p, w — w') 



(B.9) 



ImS(p,a; + W — u)') 

where the integral / dx/ [(a + + b'^\ = arctan ((a + x) /b) was used. 
In the Umit q the identity 

lim -J— e— 'VVt ^ ^s{u;') , (B.IO) 

9^0 2qVT 

allows us to perform the frequency integration, 

Im a; + iJ) = d, ^^^^ |^du; exp 

' (p + 9)^ — 1^ — M + Re ^(P; ^ 



2T; 2gVr 



arctan 



— arctan 



ImS(p, w + i^ — w') 

{p — q)"^ — Lu - fi + Re S(p, w — lu') 



Im S(p, uj + iS ~ uj') 
( {p + qY - ijJ - ^J. + Rb E(p, w - w') ' 



arctan 



\ ImE(p, w + i(5 — w') 

(p — g)^ — w — /X + Re I](p, w — uj') 



— arctan 



Im E(p, UJ + \5 — u)') 



2k^T f°^^ q 
dq 



T^P Jo [9^ + k'^? 

{p + qf -id- + Re E(p, Lo) 



arctan 



arctan, -r ^, .c-x 

lmE(p,a; + w) 

Using the following power expansion of the arctan- function 



ImE(p,w + \5) 

(p - g)2 - w - yu + Re E(p, w) 



(B.ll) 



arctan(l + a;) = - + - + 0{x'^) , (B.12) 



I.e. 

'{p + q)"^ M + s(p, uj) 



arctan 



Im E(p, UJ + iS) 

{p - q)^ - UJ - iJ + Re E(p, uj) 



arctan 



Im E(p, w + 16) 
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= 4pq |lml](p, w + i(5) 
we obtain 



- — ^ + Re i;(p, 
Imi;(p, w + iS) 



(B.13) 



Im to + i6) 



Jo 



dq 



x-ipq < Inil](p,a; + 



— a; — ^ + Re S(p, w) 
IniE(p, w + 



(B.15) 



which can be turned into 
i2 , r„2 



[Im E(p, + i(5)] + - - ^ + Re E(p, w)] 



2 Sk^T 



TT 

2kT 



dq 



(q2 + k2)2 



(B.16) 

To solve this single equation for the two unknown Re I](p, and Im I](p, uj + iS), we 
make use of the spectral representation of the Green function 

duj A{p, oj) 



doj ImE(p,a; + iS) 



, 27r z — uj 



(ImS(p,c^ + iJ))^ 

27r z — UJ 

dw ImE(p,tj + 15) 1 _T.{p,z) 



f,~p'-Re^ip,iu)y 



27r 



2kT 



2kT 



(B.17) 



(B.18) 



In the last step we also used the spectral representation of the correlated self-energy. 
The last equation has the solution 

-,1/2 



E(p,z) 



z — p- + p, 



± 



2 , N 2 
Z — p + fi 



2kT 



(B.19) 



With z = bj-\-\5 and having in mind that Im I](p, lo + iS) < for 5 > 0, we finally find 



E(p,w + i(5) 



i.e. equation (j28 



+ /i - 



sign(u; + pL-p^ 



cij + /i + iJ — 



2«;T 



1/2 



,(B.20) 



Appendix C. Details for the quasi-particle self-energy 

We start from ([25]) for the imaginary part of the self-energy using the Born 
approximation for the screened interaction potential: 



Imi;(p,cj + i&) = - 



Vtk2 

7r3/2 



dcos( 



dq 



du' ^ 

oo (<Z2 + k2)^ 



xA(p-q,tj-a;') exp exp 



,'2 



4q2T 



(C.l) 
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By assuming a frequency and momentum independent self-energy I](p,w) = 
12{0, Ef)/h), this becomes 

ImS(0,i?o/fi) = ^-^ / dcos0 / dq du' 

Since the self-energy is assumed to be independent of the frequency, the real part of the 
correlated self-energy vanishes exactly. For the Hartree-Fock part of the self-energy is 
proportional to nK^ in the classical limit [8], we also neglect this term, since it gives 
contributions of higher order in n, whereas we are only interested in the lowest order. 

After eliminating ImE(0,i?o/ft) on both sides, performing the trivial integration 
over the angle 9, which yields a factor 2, the frequency integration is performed by 
the help of [51] 

exp(-t2)di TT / 2\ r / • X o^ 
exp(-z'^) erfc(-iz) , (C.3) 



z — t 2\z 



leading to 



^..2 r dqq [°° exp(-w'Wr) dLj' 



= -A fL 2 r dgg ttRc [exp(-z2)crfc(iz)] 

V ^3'* Jo (g2+K2)2 lmE{0,Eo/h) ' ^ ■ ' 

which is rewritten as 

Im S(0, Eo/h) = -4 Y^/t' ^ {q2'^+l2)2 [exp(-^') erfc(iz)] , (C.6) 
with 

q^ + ilmJ:{0,Eo/h) 

z= . (C./j 

2g\/T 

Equation (jC.6p is given as ((55)) in section O It should be noted at this point that the 
integral converges only for finite k, i.e. the Coulomb limit k — > does not yield a 
finite result. 

Most contributions to the integral stem from small values of the wavevector, q k. 
Therefore, we may neglect the real part of z and write z = il-m'S{0, Eo/H)/2q\/T. 
Using the expansion 

lim exp(x2) erfc(a;) = - — ^ + 0{x~^) (C.8) 

in lowest order only, the q-integral can be performed, resulting in 
..2 / dqq 



exp(-(ilml](0, £;o/ft)/2(7VT)2) erfc(i • iImE(0, Eo/h)/2qVf) 



xRe 



xRe 



T 2 p dqq 

{q^ + k2)2 

exp((-Ini I](0, Eo/h)/2qVTf) erfc(-Im E(0, Eq/H) /2qVf) 
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K 



,2 



/ 

^0 



■ oc 



{q^ + «2)2 y^(-ImE(0, Eo/h)/2qVT) 
/■°° dqq^ 8Tk2 tt 



dqq 1 



(C.ll) 



(C.12) 



nlmJ:{0,Eo/h) 



(g2 + k2)2 Trim I](0, £^0/^1) 4k 



Finally, from the last line, 

lml](0,£;o/ft) 



(C.13) 



is obtained, and, after re-establishing SI units, equation ([38)) . 
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